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CHAPTER  1 
INTRODUCTION 

Among  various  kinds  of  mechanisms ,  there  is  a  certain  group  called 
force-balancing  mechanisms.  This  report  aims  to  develop  a  general  purpose 
theory  and  a  computer  code  for  design  optimization  of  force-balancing 
mechanisms  and  to  provide  example  calculations. 

1-1  Force-Balancing  Mechanisms 

A  force-balancing  mechanism  has  a  device  or  counterweight  that  main¬ 
tains  the  mechanism  in  a  force-balance  state,  or  at  least  minimizes  the 
force  necessary  to  drive  the  mechanism  over  a  certain  range  of  inputs. 
Figure  1-1  shows  an  example  of  a  force-balancing  mechanism  that  couples  an 
automobile  hood  to  the  body.  It  is  a  6  bar  linkage  with  a  linear 
mechanical  spring.  One  can  open  or  close  the  hood  easily  because  the 
spring  of  the  mechanism  exerts  a  balancing  force  against  the  weight  of  the 
hood. 

Figure  1-2  shows  a  schematic  diagram  of  an  elevating  mechanism  of  a 
weapon  with  an  equilibrator.  Preponderance  of  the  tipping  parts  is  in¬ 
evitable  in  some  applications,  hence  an  equilibrator  is  essential  for 
minimizing  elevating  or  depressing  effort.  The  weight  moment  increases 
as  the  tipping  parts  are  lowered,  so  the  equilibrator  is  designed  to  exert 
a  larger  balancing  force  when  the  elevation  angle  is  small  than  when  it  is 
large. 

Figure  1-3  shows  a  sliding  door  with  a  torsional  spring  equilibrator. 
The  counterweight  of  an  elevator  system  is  another  example  of  a  force¬ 
balancing  mechanism. 


\  m 


Torsional  Spring 


Sliding  Door 


Figure  1-3  Sliding  Door  with  a  Torsional  Spring  Equilibrator 

1-2  Static  Force-Balancing  Mechanism 
As  described  briefly  in  the  previous  paragraph,  the  force-balancing 
mechanisms  are  intended  to  maintain  parts  in  a  desired  position.  Although 
the  balancing  force  is  exerted  during  motion  of  the  mechanism,  the  motion 
is  slow.  Therefore,  the  force-balancing  mechanism  can  be  termed  a  static 
force-balancing  mechanism.  Once  the  range  of  motion  is  discretized  (see 
par.  4.2),  then  the  mechanism  is  analyzed  at  each  input  point  as  a  static 
case.  Hence,  the  term  static  force-balancing  mechanism  is  quite  natural. 

1-3  Linear  Equilibrator 

To  avoid  being  too  broad,  this  report  deals  only  with  static  force¬ 
balancing  mechanisms  that  have  linear  spring  equilibrators.  As  the  word 
"linear"  implies,  the  equilibrator  has  linear  characteristics.  The 


equilibrator  may  be  a  linear  mechanical  spring,  pneumatic  spring,  or 
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torsional  spring.  The  magnitude  of  the  force  exerted  by  the  equilibrator 
is  proportional  to  the  change  of  displacement;  i.e.,  the  balancing  force  is 


F  =  K  AS,  (1-1) 

where  F  =  balancing  force  or  torque 
K  =  constant  spring  rate 
A£  =  linear  or  angular  displacement 
In  most  cases,  the  force  necessary  to  equilibrate  a  mechanism  is  non¬ 
linear,  while  the  equilibrating  force  available  is  linear.  Therefore, 
perfect  balancing  is  impossible  over  the  entire  range  of  motion.  This  fact 
leads  to  a  problem  of  worst  case  optimization  (see  Ref.  1,  Ch.  3,  or  Ref.  2), 
which  is  closely  related  to  min-max  theory  (see  Ref.  1). 

To  solve  the  equations  for  determining  the  state  of  a  mechanism,  a 
technique  based  on  the  DADS  formulation  (see  Ref.  2)  is  used.  Using  the 
DADS  formulation,  one  can  minimize  the  computational  difficulty  of  solving 
highly  nonlinear  equations. 
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CHAPTER  2 

STATIC  ANALYSIS  OF  MECHANISMS 

For  obtaining  the  equations  determining  the  state  of  a  mechanism,  a 
constrained  multi-element  formulation  is  described  in  this  section.  Kine¬ 
matic  constraint  equations  for  the  most  common  joints  used  in  mechanisms 
are  derived  as  state  equations. 


2-1  Constrained  Multi-Element  Formulation 
Figure  2-1  shows  the  ith  body  of  a  mechanical  system  that  lies  in  an 
inertial  x-y  reference  frame.  A  local  coordinate  system  n^)  is  em¬ 

bedded  in  the  body.  The  position  of  the  local  coordinate  system  is 
arbitrary,  but  is  usually  located  at  the  center  of  mass.  The  position  and 
orientation  of  the  body  are  described  by  three  generalized  coordinates; 
x±,  y±,  and  0i. 

The  constrained  multi-element  formulation  regards  joints  connecting 
two  bodies  as  constraints  imposed  on  the  relative  motion  of  the  connected 
bodies.  By  doing  this,  the  method  is  capable  of  modeling  open  or  closed 
loop  systems.  However,  because  each  joint  (revolute  joint  and  translational 
joints  are  most  common)  brings  forth  two  additional  constraint  equations,  it 
gives  rise  to  a  larger  number  of  equations  than  other  methods. 


2-2  State  Equation 

Consider  a  system  with  n  bodies  connected  by  Z  joints  in  the  plane. 

Each  body  has  three  generalized  coordinates,  x,  y,  and  0.  Hence  a  system 
of  n  bodies  has  a  total  of  3n  generalized  coordinates.  As . stated  in  the 
previous  paragraph,  a  joint  is  regarded  as  a  constraint  on  the  motion  of  the 
pair  of  bodies  that  the  joint  connects.  The  most  common  joints  (revolute 
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and  translational)  are  dealt  with,  each  joint  having  two  constraint  equations. 
The  constraint  equations  express  necessary  conditions  that  the  two  bodies 
must  satisfy  to  be  compatible  with  the  joints.  A  system  that  has  5,  indepen¬ 
dent  joints  has  2 5,  constraint  equations.  Thus,  the  number  of  free  degrees- 
of-freedom  of  this  system  is  m  =  3n  -  25,.  The  condition  m  >  0  must  be 
satisfied  by  all  kinematic  systems.  The  3n  degrees-of-freedom  are  denoted 
by  z  e  R  ,  z  =  [x^,  y^,  9^,  ...»  x^,  y  ,  0  ]  .  Since  the  components 

of  the  vector  z  describe  the  state  of  the  system,  they  are  called  state 
variables  and  z  is  called  state  vector.  The  25,  constraint  equations  are 
the  kinematic  conditions  that  the  state  must  satisfy,  thus  these  constraint 
equations  are  called  kinematic  constraint  equations  and  can  be  written  as 

$k(z,  b)  =  0  (2-1) 

where 


$k(z,  b) 


<|>k(z,  b) 
^25,(z ’  b) 


and  b  e  R  is  a  vector  of  design  variables. 

The  basic  purpose  of  designing  a  kinematic  system  is  to  obtain  a  system 
that  allows  the  transmission  of  motion  or  force  from  an  input  link  to  an 
output  link.  The  mechanism  can  only  be  given  input  motion  through  a  set  of 
free  degrees-of-freedom,  which  is  specified  as  a  function  of  some  free 
parameter,  or  as  some  relationship  between  the  3n  state  variables.  Since 
all  degrees-of-freedom  must  be  specified  to  drive  the  mechanism,  m  additional 
driving  equations  of  constraint  are  required,  in  the  form 


$^(z,  b,  a)  =  0 


(2-2) 
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where 


4>d(z,  b,  a) 


4>d(z.,  b,  a) 
<J>d(z,  b,  a) 


and  a  e  is  a  vector  of  environmental  parameters  (see  Ref.  1,  Ch.  3). 

The  range  of  motion  is  specified  by  a.  For  example,  consider  a  four-bar- 
linkage  of  Fig.  2-2.  This  mechanism  has  one  free  degree-of-freedom. 
Hence,  a  could  be  the  angle  of  the  input  link  with  respect  to  the  ground. 
Thus  a  specifies  the  position  at  which  the  mechanism  is  being  analyzed. 


W77j777777777/777777777777777777177W7Z 


Figure  2-2  Four-Bar-Linkage  with  One  Degree-of-Freedom 


Combining  Eqs.  2-1  and  2-2,  one  has  3n  constraint  equations,  which  can 
be  written  in  the  form 


4>(z,  b,  a) 


4>k(z,  b) 
$d(z,  b,  a) 


(2-3) 
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Equation  2-3  is  the  state  equation  for  position  analysis  of  the  mechanism. 
Of  all  the  variables  in  Eq.  2-3,  design  variables  b  and  environmental  para¬ 
meters  a,  are  specified  before  any  attempt  is  made  to  solve  the  equations 
for  the  state  z.  Thus,  Eqs.  2-3  are  interpreted  as  3n  equations  with  3n 

unknowns  z.  This  system  can  be  solved  uniquely  for  z(b,  a),  provided  that 
the  implicit  function  theorem  (Ref.  3)  is  satisfied;  i.e.,  the  matrix 


2i 

3z 


3nx3n 


is  nonsingular. 


2-3  Kinematic  Constraint  Equations 
for  Typical  Joints 


2-3.1  Revolute  Joint 

Figure  2-3  shows  two  adjacent  bodies  (i  and  j)  that  are  connected  by 

revolute  joints.  The  local  coordinate  systems  (£  ,  n.)  and  (£  ,  n.)  are 

ii  j  J 

embedded  in  each  body  and  located  by  x^,  y^,  0^  and  ,  y^ ,  0  ^ ,  with  respect 

to  the  global  (x,  y)  coordinate  system.  The  origins  of  the  local  coor- 

dinates  are  positioned  by  R^  and  R ^ ,  with  respect  to  global  coordinates. 

Consider  two  points  P  .  on  the  ith  body  and  P  .  on  the  jth  body  that 

ij  3 1 

are  located  by  r^  and  r  ^  with  respect  to  the  local  coordinate  system. 

The  connection  of  P^  and  P^  is  expressed  by  a  vector  r^.  Executing  a 
closed  loop  from  the  origin  of  the  global  coordinate  yields  the  vector 
relationship 


Ri +  rij  +  rP  -  rji "  Ej ' 0 


(2-4) 


The  kinematic  constraint  equations  for  a  revolute  joint  can  be  obtained  by 

requiring  that  r"  =0;  i.e.,  requiring  that  P.  coincide  with  P...  Using 

P  ij  Ji 
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the  proper  coordinate  transformation,  Eq.  2-4  can  be  written  as 


<t>x  »  Xi  +  5^cos 0^  -  n^ine^^  -  Xj  -  ^cosS 

+  rijSin6j  =  0 

<{>y  =  Y.^  +  f^sinG^^  +  r^cosG^.  -  -  ^sinG^ 


-  n j  cos  0^  «*  0 


(2-5) 


2-3.2  Translational  Joint 

Figure  2-4  shows  the  geometry  of  a  translational  joint.  Definitions 
of  the  local  coordinates  and  position  vectors  are  similar  to  those  of  a 
revolute  joint  in  the  previous  paragraph.  The  translational  joint  is 
defined  by  requiring  that  the  vector  r^  be  parallel  to  the  path  of  relative 
motion  between  the  two  bodies ,  where  r^  is  the  vector  between  two  points 

.  ^  y 

P . .  and  P,.  that  are  located  by  two  nonzero  vectors  r . .  and  r.,  that  are 

ij  ji  ij  ji 

perpendicular  to  the  line  of  relative  motion.  From  Eq.  2-4, 


v-<^+ vw 


_y. 

The  condition  that  the  vector  r„  is  perpendicular  to  the  vector  r^  gives 
one  constraint  equation  as 


r  •  r  =  0 
P  ij 


(2-6) 


-y. 

The  second  constraint  equation  is  obtained  by  noting  that  r, ,  and  r..  are 

ij  3i 

parallel  to  each  other;  i.e., 


riJ  x  rji  ’  * 


(2-7) 


Equations  2-6  and  2-7  can  be  rewritten  in  scalar  form  as 
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(ui  -  X1)(D1  -  v  +  <V1  -  V<vi  -  y  - 0 
+9  -  <Ui  -  Xi><Vj  -  V  -  »±  -  V(UJ  '  V  ’  0 

where 

+  £^cos0^  -  ri^sin©^ 

Uj  =  Xj  +  SjCoaOj  -  rij3in0j 

Vi  =  +  S^sinO^  +  nicos0i 

V.  =  Y.  +  £.sin0  +  n. cos  0. 
j  j  j  J  3  j 


(2-8) 


Equations  2-8  are  the  constraint  equations  for  the  translational  joint. 
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CHAPTER  3 

DERIVATION  OF  GENERALIZED 
FORCE  FOR  A  LINEAR  SYSTEM 


Consider  a  system  of  n  bodies  that  has  m  linear  equilibrators  and  s 

torsional  spring  equilibrators.  A  system  that  has  linear  force-balancing 

elements  is  called  a  linear  system  in  this  Chapter.  Figure  3-1  shows 

adjacent  ith  and  jth  bodies,  with  the  kth  linear  equilibrator  element  and 

the  qth  torsional  equilibrator  element.  The  definition  and  location  of 

local  coordinates  of  each  body  are  similar  to  those  for  a  revolute  joint. 

However,  in  this  system  the  kth  linear  spring  equilibrator  is  connecting 

the  ith  and  jth  bodies.  The  attachment  points  and  P  ^  are  located  by 

vectors  r^  and  r^*  with  respect  to  each  body  fixed  coordinate  system. 

The  qth  torsional  spring  equilibrator  is  installed  in  ith  body.  The  length 

of  the  linear  spring  equilibrator  is  expressed  as  £  .  The  linear  spring 

I  jk. 

-fc- 

equilibrator  forces  are  expressed  as  and  F  The  torque  exerted  by  a 

torsional  spring  equilibrator  is  expressed  as  M_^. 

Similar  to  Eq.  2-4,  executing  a  closed  loop  from  the  origin  of  the 
global  coordinate  system  yields  the  vector  relationship 


4-  ~y  -y 

R.  +  r.,  +  £..,  -  r.,  -  R  =  0 
i  ik  ijk  jk  j 


(3-1) 


From  Eq.  3-1,  £„k  can  be  written  as 


->■  *4-  -V 

£...  =  -R.  -  r. ,  +  R.  +  r.. 
ijk  i  ik  j  jk 


(3-2) 


The  derivation  of  generalized  force  starts  from  the  principle  of 
virtual  work  (see  Ref.  4);  i.e. , 


6W  =  0 


(3-3) 
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The  virtual  work  due  to  the  applied  force  during  the  virtual  displacement  is 
6W  =  FT  fir  (3-4) 


where 


F  =  applied  force 
fir  =  virtual  displacement 

For  the  force  F,  the  virtual  displacement  of  the  point  of  application  of  the 
force  is 


.  9<Ei +  V  . 

6r - jj -  Sz 


(3-5) 


Let  the  virtual  work  done  by  equilibrators  and  by  external  forces  be 

5W  and  6W  ,  respectively.  Then  the  virtual  work  done  by  r  linear 
eq  ex 

spring  equilibrators  and  s  torsional  spring  equilibrators  is 


*  l  9(R  +  r  ) 

6W  =  l  \  F 1  - ~z — —  fiz  +  F, 


eq  i  ik 


3z 


T  +  r1k> 

jk  3z 


s  36. 

+  y  M.  r-=>  fiz 
q-l  92 


(3-6) 


The  virtual  work  done  by  external  forces  and  torques  is 


fiW 


ex 


n 

-  I 

i=l 


g  3(R  +  r  )  h  30 

?  P* - fiz  +  Y  T.  -^-5  = 

„ i£  3z  L.  nn  3z 

1=1  m=l 


(3-7) 


where 

fiW  =  virtual  work  by  external  forces 
ex 

P.  =  £th  external  force  on  body  i 
i£ 

T.  =  mth  external  torque  on  body  i 
ira 

g  =  number  of  external  forces  on  body  i 

h  =  number  of  external  torques  on  body  i. 


The  total  virtual  work  is  thus 


SW  =*  6W  +  6W 

eq  ex 


(3-8) 


Farkas  Lemma  (see  Ref.  1)  guarantees  the  existence  of  a  Lagrange 
multiplier  vector  y  e  R^n  such  that 


6W  -  y  6$  =  0 


(3-9) 


where 


$  =  the  vector  of  state  equations. 


Substituting  6W  and  6Wgx  in  Eqs.  3-6  and  3-7  in  Eq;  3-8  and  then  6W  in 
Eq.  3-8  into  Eq.  3-9  and  factorizing  <Sz, 


f  iF  T  3(R1  +  rlk)  f  F  T  ,3(Rj  Uikj 
^  |  ik  3z  jk  3z 


A  (  g  3(R.  +  r  )  h  30 

+  ix  J/- 


s  38.  _ 

q-i  J 


(3-10) 


Since  Eq.  3-10  holds  for  arbitrary  §z. 


?  i  3(Ri  +  ra)T  r  ,  atRi  +  rik)T  r 

l  )  az  tk  az 


n 

+  I 

i=l 


g  3  (R .  +  r  )T  h  30  T 

y  — i - p  +  y  _A_ 

L  az  c  1  az 


iH  +  ^  32  Tim 

m=l 


s  36, 

+  I,-jrMiq 

q=l 


Tz“y 


(3-11) 


Defining  the  left  hand  side  of  Eq.  3-11  as  a  vector  G,  Eq.  3-11  can  be 


written  as 
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8$  (z,  b,  a) 

- - -  p  =  G(z,  a,  b) 


(3-12) 


where 


r  (3(R.  +  r,,  )  8  (R,  +  r,,  ) 

G(z,  a,  b)  =  l  l - i- - — —  F . .  +  - I- - li-F 

*  k=i  )  92  ik  82  3k| 


s  80. 

+  l  — r^—  M . 

q=l  8z  ^ 


n 


+  I 


3<Ri+ri4>' 


i=l  £=1 


8z  14. 


h  80. 

p-+  l  (3-13> 

m=l 


By  Newton's  law  of  motion,  =  -  F_.^.  Hence  substituting  F.^  into 
Eq.  3-13,  and  denoting  the  sum  of  the  first  two  terms  of  Eq.  3-13  as 
Geq(z,  b,  a), 


Geq(z,  b,  a)  =  £ 


r  ^(Ri  +  rik)T„  9(Rj  +  rik)T 


k=l 


8z  Fik  -  -  8z  J  Fik! 


s  80. 

+  I  — r~—  M, 
q=l  9z  iq 


(3-14) 


Factoring  F^>  this  is 


G  (z,  b,  a)  =  l 


r  |  +  rlk)T  3(R.  +  rik)T 


eq 


k=l 


8z 


8z 


ik 


T 

s  80. 

+  I  M. 

q=l  32  1£1 


(3-15) 


which  becomes 
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Ceq(z,  b,  a) 


where 

G  =  rigidity  modulus 
d  =  wire  diameter 
D  =  coil  diameter 
n  =  nunfcer  of  effective  coils 

In  a  linear  mechanical  spring  equilibrator,  a  compression  spring  is 
generally  used,  rather  than  a  tension  spring,  because  no  extra  ending  is 
necessary  for  installation  and  if  failure  of  the  spring  occurs  the  spring 
will  still  sustain  a  reduced  load.  Figure  3-2  shows  a  pull-type  equilibrator, 
which  has  a  compression  spring.  Figure  3-2 (a)  shows  an  equilibrator  spring 


*.q  =  spring  free  length 
i  B  cylinder  length 
ip  =  piston  length 
i  =  equilibrator  length 


Figure  3-2  Spring  Deflection  Calculation 
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before  assembly.  The  spring  free  length  is  usually  longer  than  the  cylinder 
length,  for  presetting.  Figure  3-2 (b)  shows  the  compressed  spring  when  the 
equilibrator  length  is  Z.  Spring  deflection  is  the  sum  of  lengths  A  and  B, 
which  are  -  Z  ^  and  Z  -  respectively;  i.e., 


AZ 


(Z0  -  zc)  +  (Z  -  y 


z  +  z 

p 


(3-19) 


where 

Zq  =  spring  free  length 
Z =  cylinder  length 
«  piston  length 
Z  =  equilibrator  length 

Equations  1-1,  3-18,  and  3-19  determine  the  magnitude  of  the  equilibrator 
force.  The  direction  of  the  equilibrator  force  is 


n  =  (3-20) 

where  n  =  direction  vector 

Exploiting  Eqs.  1-1,  3-18,  3-19  and  3-20,  the  first  term  on  the  right 
side  of  Eq.  3-17  becomes 


where  k  represents  the  kth  equilibrator. 
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CHAPTER  4 

OPTIMAL  DESIGN  PROBLEM  FORMULATION 


4-1  Statement  of  the  Optimal  Design  Problem 

The  optimal  design  problem  for  a  force-balancing  mechanism  can  be  stated 

as  follows : 

k. 

Find  b  e  R  to  minimize  the  cost  function 

*0  =  *0(b>  (4-1) 


subject  to  the  constraints; 

(i)  state  equation  for  position,  from  Eq.  2-3, 
$(z,  b,  a)  =  0 

(ii)  force  equation,  from  Eq.  3-12 

T 

9$  (z.  b ,  a)  „  , 

-  *-*=-* — "  V  =  G(z,  a,  b) 

(iii)  composite  design  constraints 

max  \p^(z,  y(b,  a))  ^  0,  i  =  1,  p 

a 

where  p  is  the  total  number  of  design  constraints. 


(4-2) 


(4-3) 


(4-4) 


4-2  Discretized  Problem 

The  optimal  design  problem  stated  in  the  previous  paragraph  is  valid  on 
the  entire  range  of  a.  Such  problems  are  called  parametric  optimal  design 
problems.  Techniques  for  handling  parametric  optimal  design  problems  have 
been  developed  in  Ref.  1.  One  simple  way  of  handling  this  problem  is  to 
discretize  a  into  finite  number  of  grid  points.  The  composite  design  con¬ 
straints  are  required  to  hold  at  every  grid  point  on  the  ot-grid,  so  the 
generalized  force  must  be  calculated  at  every  grid  point.  Thus  the  state 
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equation,  force  equation,  and  composite  design  constraints  can  be  rewritten 
as 


$(a,  b,  orb  =  0, 

j  —  i»  •••»  t 

(4-5) 

3$(z,  b,  qJ)T  _ 
dz  p 

G(z,  b,  orb  ,  j  =  1 . t' 

(4-6) 

max  $  (z,  y(b ,  orb,  b)  <  0,  i  =  1 . .  (4-7) 

j 

j  =  1»  •  •  •  i  t* 

where  t'  is  the  total  number  of  grid  points  of  a. 

4-3  A  Reformulation  of  the  Problem 

It  is  now  necessary  to  define  the  cost  function  more  precisely.  In  a 
static  force-balancing  mechanism,  the  eventual  aim  is  to  minimize  the 
deviation  between  the  balancing  force  and  the  force  to  be  balanced.  Since 
this  imbalance  varies  as  the  environmental  parameter  a  varies ,  the  maximum 
unbalance  is  chosen  to  be  minimized.  One  can  now  introduce  an  artificial 
design  parameter  b^+^  to  satisfy  the  constraint 


max  f(z,  p(cr*,  b)  ,  b) 

j 


'  bk+l  *  ° 


j  -  1, 


t' 


(4-8) 


where  f  is  the  unbalance  between  the  balancing  force  and  the  force  to  be 
balanced.  This  upper  bound  on  the  maximum  unbalance  can  be  modified  as 


Vb) -  bk+i  (4"9) 

with  the  additional  constraint  of  Eq.  4-8. 

Since  the  range  of  a  is  being  discretized,  the  constraints  must  be 
evaluated  at  each  a-grid  point  to  determine  violations  in  the  constraints. 
The  relative  maximum  strategy  (see  Par.  4-5)  is  used  for  the  constraint 
equations  of  Eqs.  4-7  and  4-8. 
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The  optimal  design  problem  may  now  be  restated  as: 
k+1 

Find  b  e  R  to  minimize  the  cost  function 


^0  bk+l 


(4-10) 


subject  to  the  constraints; 

(i)  state  equation  for  position,  from  Eq.  4-5, 


$(z,  b,  or*)  =  0,  j  =  1,  ...,  t’ 


(4-11) 


(ii) 


force  equation,  from  Eq.  4-6, 


9$(z,  b,  aV 
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y  =  G(z,  b ,  or ) ,  j  =  1,  . . .  ,  t' 


(4-12) 


(iii)  composite  design  constraint 


<K(z,  y (b ,  or3),  b)  <  0,  i  =  1,  p 


(4-13) 


j  -  1 ,  •  •  •  ,  t ' 


where  t'  is  the  total  number  of  grid  points  of  a. 


4-4  Design  Sensitivity  Analysis 

Derivatives  of  the  cost  function  and  constraint  equations  with  respect 
to  design  variables  must  be  calculated  before  any  optimization  algorithm  is 
applied.  The  idea  of  design  sensitivity  analysis  is  to  construct  an 
approximation  of  the  design  problem  that  can  be  analyzed  to  determine  the 
effect  of  a  change  6b  in  b^.  The  approximate  problem  is  obtained  by  making 
linear  approximations  to  nonlinear  functions.  A  linear  approximation  to 
the  change  in  is  trivial  and  linear  approximations  to  changes  in  ip^,  due 
to  the  small  changes  in  all  variables ,  are 

9<K  dip 

6+i  =  IT  62  +  —  6y  +  InT 


6b 


(4-14) 
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To  construct  a  design  sensitivity  vector  is  to  find  the  vector  as 


5^i  =  6b 


(4-15) 


Hence,  it  is  necessary  to  express  Sz  and  6p  in  terms  of  6b.  To  get  the 
relation  between  6z  and  6b,  take  the  first  variation  of  Eq.  4-11,  with  a 
fixed,  to  obtain 


3?  Sz  +  3b  5b  "  ° 


(4-16) 


Since  —  is  non-singular  (see  Par.  2-2), 

d  Z 


-  (-  f)'1  f  «• 


(4-17) 


The  first  term  of  Eq.  4-14  thus  becomes 


3^  62  “  "  li 


i  /_a#y 

\3z/ 


(4-18) 


Defining  an  adjoint  variable  X^  such  that 


X!T  =  _  / 3*V 

l  3z  \3z/ 


(4-19) 


reduces  the  cost  of  computation  by  avoiding  calculation  of  the  inverse  matrix 
/3$\-^  1 

;  i.e.,  X^  is  the  solution  of  adjoint  equation 


a if  a  .  (!!iV 

\3z /  i  \3z  / 


(4-20) 


Now,  Eq.  4-18  becomes 


_  1  3$ 

—  6z  -  Xt  w  6b 


(4-21) 


26 


One  can  also  express  6y  in  terms  of  6b  by  taking  the  first  variation  of 
Eq.  4-12,  keeping  b  invariant. 


\3z  / 


x  3G  .  6G 

611  '  ta  62  +  6b  Sb 


(4-22) 


This  equation  yields 


-  '0\' 


[(If)1*] 


6z  +  6b 
ob 


The  second  term  of  Eq.  4-14  thus  becomes 


3y  6y  3y 


i  a  [(if)1-] 

l  \  3z  /  J  )  3z  3z 


6z  +  ~  6b 
ob 


Similarly  to  Eq.  4-19,  define  an  adjoint  variable  xf  as 


(4-23) 


(4-24) 


■  ©  f@)' 


(4-25) 


This  can  be  rewritten  as 


ii  x2 .  (S) 

3z  Ai  \3y  / 


(4-26) 


Using  Eq.  4-17  and  4-25,  Eq.  4-24  becomes 


9*i  2T 

—  6y  "  Xi 


[mT»] 


/_3i\-1  ii 

\  3z/  3b 


6b  +  |~  6b 
9b 


(4-27) 


Defining  such  that 


[(If)1-] 


m' 


(4-28) 
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one  can  rewrite  Eq.  4-27  as 


3*i  2T  3T  9$  9G 

9T  6y  =  Xi  Xi  lb  +  3b  5b 


and  X.  is  the  solution  of 
x 


(If) 


[(If) 

T  ' 

l 

9z  ( 

Substituting  Eqs.  4-21  and  4-29  into  Eq.  4-14,  one  finally  has 


if  ,1  '  0T  ,T“|  0T  ^ 

.1  i.  ,  2  .3  9$  ,  ,2  9G  , 

X.  +  X.  X.  —  +  X.  —  + 

Li  i  x  J  9b  19b 


6*i  -  ^Xj  +  X~  X.  J  - 
The  design  sensitivity  vector  is  thus 


‘iT=  [xi +  xf  xf]  (ff)+x 


2X  j>G  ^ i 

9b  +  9b 


(4-29) 


(4-30) 


8*i  )  iT 

~~  >  6b  =  6b  (4-31) 


(4-32) 


12  3 

where  the  adjoint  variables  X^,  X^,  and  X^  are  the  solutions  of  Eqs.  4-20, 
4-26,  and  4-30,  respectively. 


4-5  Optimization  Algorithm 

With  the  design  sensitivity  vector  calculated  in  the  previous  paragraph, 
one  can  implement  an  optimization  algorithm  of  his  choice.  The  gradient 
projection  method  with  constant  error  correction  (Ref.  1)  is  chosen  here. 

Among  the  constraint  equations  of  Eq.  4-13,  some  are  dependent  on  the 
environmental  parameter  a.  Such  constraints  are  called  parametric  con¬ 
straints.  Constraints  that  are  independent  of  the  environmental  parameter 
are  called  nonparametric  constraints.  This  distinction  is  necessary  to 
save  calculations,  because  one  does  not  need  to  calculate  nonparametric 
constraints  at  every  grid  point.  Enforcing  constraints  only  at  relative 
maxima  in  a,  called  the  relative  maximum  strategy,  decreases  the  number  of 
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calculations  even  further.  This  strategy  chooses  the  grid  point  at  which 
the  relative  maximum  violation  of  a  constraint  occurs  and  solves  the  state 
equations  only  at  that  point  and  in  a  neighborhood  of  the  point.  This  is 
called  a  non-sweep  iteration.  The  actual  computational  optimization 
algorithm  is  presented  in  Chapter  3  of  Ref.  2  as  follows: 

Step  1  :  Estimate  a  design  and  impose  a  grid  on  the  range  of  the  input 
parameter  a. 

Step  2  :  Solve  the  state  equations  of  Eq.  4-11  and  4-12  for  z(ci'')  and 
hOr1),  respectively. 

Step  3  :  Check  all  non-parametric  and  parametric  design  constraints.  Form  a 

reduced  vector  \p  consisting  of  all  e-active  non-parametric  constraints 

and  relative  maxima  of  e-active  parametric  constraints. 

12  3 

Step  4  :  Compute  the  adjoint  variables  ,  A„  ,  and  A^.  from  Eqs.  4-20, 

4-26,  and  4-30,  associated  only  with  the  constraint  functions  in  ij;. 
iT 

Construct  £  of  Eq.  4-32  for  each  of  these  constraint  functions. 

i’  iT 

Form  the  matrix  A  =  [£  J],  whose  columns  are  the  vectors  SL 

~  ~  T 

corresponding  to  constraint  functions  in  ip.  Thus,  Sip  =  A  6b. 

Step  5  :  Compute  the  vector  M  and  the  matrix  from  the  following 
relations : 


M  =  AT  W"1  1° 

**o 


(4-33) 


where  sP  is  obtained  from  the  variation  of  cost  function  as 


8*0  0T 
%  “  IT  6b  B  z  6b 


(4-34) 


T  -1 

M,  =  Aa  W  A 


(4-35) 
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CHAPTER  5 
APPLICATIONS 

5-1  A  Linear  Mechanical  Spring  Weapon  Equilibrator 
As  an  example  calculation,  a  linear  mechanical  spring  equilibrator  for 
a  weapon  is  designed.  Two  concentric  springs  are  installed.  Their 
direction  of  winding  is  opposite,  to  prevent  the  inner  spring  from  inter¬ 
lacing  with  the  outer  spring.  Both  cylinder  wall  and  piston  rod  work  as 
guides  to  prevent  buckling. 

Figure  5-1 (a)  shows  geometry  of  the  equilibrator  and  the  definitions  of 
global  and  local  coordinates.  Figure  5-1  (b)  shows  dimensions  in  the  initial 
position.  Body  I  (top  carriage)  is  the  weapon  support  and  body  II  (tipping 
parts)  is  connected  to  body  I  at  the  trunnion  (a  revolute  joint).  This 
system  has  only  one  moving  body,  hence  the  total  number  of  degrees-of-freedom 
is  3.  One  revolute  joint  gives  rise  to  two  state  equations  in  Eq.  2-5. 
Therefore  this  system  has  one  free  degree-of-f reedom  and  has  one  driving 
constraint,  which  is  written  in  the  form 

0  -  dj  =  0  j  =  1,  t'  (5-1) 

where 

0  =  elevation  angle 

ou  =  given  values  at  grid  point  j 

t*  =  total  number  of  grids 

Equations  2-5  and  5-1  form  the  state  equations  (Eq.  4-11). 

The  inner  and  outer  springs  are  installed  in  the  same  equilibrator, 
hence  there  is  only  one  equilibrator  length  vector  £.  The  force  equation 
of  Eq.  3-21  becomes 


Top  Carriage  (I) 


I 

L 

(a)  Geometry  of  Equilibrator 


_ \ 


(-0.1543,  -0.2673) 


(b)  Initial  Condition  of  Equilibrator 

Force  in  N,  Length  in  m  and  Angle  in  Radian 


Figure  5-1  Geometry  of  Weapon  Equilibrator 
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T 

=  jHL 

az 


Gd, 


8Dl"! 


£  -  £  + 
c  p 


+ 


Gd, 


8D2n2 


(5-2) 


The  design  variables  are  chosen  as  spring  wire  diameters,  numbers  of 
effective  coils,  and  pitches;  i. e. 


b..  =  d. 


b„  =  n 


b  = 


b,  = 


1 

Pi 

d„ 


bc  =  n 


b,  = 


2 

P2 


(5-3) 


where  d^,  d 2  =  wire  diameters  of  spring  1,  2 

nl*  n2  =  num^er  °f  effective  coils  of  springs  1,  2 
pl’  P2  =  Pitches  sPring  2 

Noting  that  the  spring  free  length  is  the  number  of  active  coils  times 
pitch  and  using  Eq.  5-3,  Eq.  5-2  becomes 


9z 


(*> 


h4 

I  12  bl 

4.18  x  10  ~  (b-  b0  +  H.  -  1.41) 
I  b2  J  Z  J 


13  b4 

+  1.41  x  10  -t2  (b ,  bc  +  Jt.  -  1.41) 

b^  6  5  j 


(5-4) 
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where  G  =  7.93  *  10^  Pa 

=  0.13335  m 

D2  =  0.889  m 

£  =  0. 7  m 

c 

£  =  0.71  m 
P 

The  inner  spring  is  required  to  support  40%  of  the  outer  spring  load. 
This  force  distribution  gives  the  constraint 


.4  h4 

b  -I  b  / 

0.1185  ~  (b3  b2  +  H.  -  1.41)  -  —  (b&  b5  +  SL.  -  1.41)  <  0  (5-5) 


where  numberical  values  of  G,  ,  D2>  £ c>  and  £p  are  substituted. 
The  shear  stress  of  a  spring  under  load  P  is  (Ref.  5) 


8PD 


g 

Allowable  stress  is  taken  as  x  =  9.1  x  10  Pa  (Ref.  5)  and 

a 

Pj^  =  F/1.4 
P2  =  0.4F/1.4 


(5-6) 


(5-7) 


are  the  loads  for  springs  1  and  2,  where 


F  =  equilibrator  force 

g 

Substituting  P,  and  P-  and  using  x  =  9.1  *  10  Pa,  the  two  stress  constraints 
12  3- 

are 

bl 

1559.8  (b,  b0  +  £,  -  1.41)  -  1  <0 

b2  3  2  J 

b , 

3508.8  jp  (b6  b5  +  £j  -  1.41)  -  1  *  0 


(5-8) 
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This  optimization  problem  may  now  be  stated  as  follows: 
Find  be  to  minimize 


* 


0 


b 


7 


subject  to  state  equations  of  Eq.  4-11, force  equation  of  Eq.  5-4,  constraint 
equations  of  Eqs.  5-5  and  5-8,  and 


(5-9) 


where  p,  is  the  constraint  force  in  the  z  direction;  i.e.,  the  direction  of 
o  o 

rotation  of  the  tipping  parts.  This  is  interpreted  as  the  unbalance  between 
the  weight  moment  of  the  tipping  parts  and  the  equilibrator  torque. 

The  starting  values  of  design  parameters  were  taken  as  shown  in  Table  5-1. 

■I  o 

This  resulted  in  |  1 6b  ||  =  160.7,  |  |6b^|  |  =  13.14,  and  a  cost  function  value 
of  51.24.  The  constraint  of  Eq.  5-9  was  violated  and  other  constraints  were 
active  at  some  grid  points.  After  11  iterations,  the  constraint  of  Eq.  5-9 

became  active.  After  20  iterations,  j|5b^(|  =  2.55  x  10 

O  __  O 

| | 6b ^ | |  =  1.46  x  10  ,  and  the  cost  function  value  was  94.29.  All  the  con¬ 

straints  are  satisfied,  but  are  active  at  some  grid  points.  The  total  com¬ 
puting  time  for  20  iterations  on  an  IBM  370/168  was  35.3  sec.  The  cost 
function  history  is  shown  in  Fig.  5-2. 


Table  5-1  Result  for  the  Weapon  Equilibrator  Design 


m 


m 


d2  m 


P2  m 


max  moment 
unbalance 
Nm 


Starting 

Values 

Optimum 

Design 

0.0165 

0.01663 

13.9 

13.84 

0.0631 

0.06247 

0.0108 

0.0109 

19.7 

19.83 

0.0426 

0.04232 

51.2* ** 

94.29 

Conventional 

Design* 


0.0171 


0.0101 


*See  Ref.  6. 

**Infeasible  design 
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US  Army  Research  Office  (02) 

P.0.  Box  12211 
ATTN:  Dr.  F.  Schmiedeshoff 

Dr.  R.  Singleton 
Research  Triangle  Park,  NC  27709 

HQ,  DA  (01) 

ATTN :  DAMA-AR 

Dr.  Herschner 
Washington,  D.C.  20310 

HQ,  DA  (01) 

Office  of  Dep  Chief  of  Staff 
for  Rsch,Dev  &  Acquisition 
ATTN :  DAMA-AR 

Dr.  Charles  Church 
Washington,  D.C.  20310 

HQ,  DARCOM 

5001  Eisenhower  Ave . 

ATTN:  DRCDE 

Dr.  R.L.  Haley 
Alexandria,  VA  22333 


Michigan  48090,  of  corrections 


Director  (01) 

Defense  Advanced  Research 
Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

Commander  (01) 

US  Army  Combined  Arms  Combat 
Developments  Activity 
ATTN :  ATCA-CCC— S 
Fort  Leavenworth,  KA  66027 

Commander  (01) 

US  Army  Mobility  Equipment 
Research  and  Development  Command 
ATTN:  DRDME-RT 

Fort  Belvoir ,  VA  22060 

Director  (02) 

US  Army  Corps  of  Engineers 
Waterways  Experiment  Station 
P.0.  Box  631 
Vicksburg,  MS  39180 

Commander  (01) 

US  Army  Materials  and  Mechanics 
Research  Center 
ATTN:  Mr.  Adachi 

Watertown,  MA  02172 

Director  (03) 

US  Army  Corps  of  Engineers 
Waterways  Experiment  Station 
P.0.  Box  631 
ATTN:  Mr.  Nuttall 

Vicksburg,  MS  39180 

Director  (04) 

US  Army  Cold  Regions  Research 
&  Engineering  Lab 
P.0.  Box  282 

ATTN:  Dr.  Freitag,  Dr.  W.  Harrison 

Dr.  Liston,  Library 
Hanover,  NH  03755 


(01) 


Grumman  Aerospace  Corp  (02) 
South  Oyster  Bay  Road 
ATTN:  Dr.  L.  Karafiath 
Mr .  F .  Markow 
M/S  A08/35 
Bethpage,  NY  11714 

Dr.  Bruce  Liljedahl  (01) 

Agricultural  Engineering  Dept 
Purdue  University 
Lafayette,  IN  46207 

Mr.  H.C .  Hodges  (01) 

Nevada  Automotive  Test  Center 
Box  234 

Carson  City,  NV  89701 

Mr.  R.S.  Wismer  (01) 

Deere  &  Company 
Engineering  Research 
3300  River  Drive 
Moline,  IL  61265 

Oregon  State  University  (01) 
Library 

Corvallis,  OR  97331 

Southwest  Research  Inst  (01) 

8500  Culebra  Road 

San  Antonia,  TX  78228 

FMC  Corporation  (01) 

Technical  Library 
P.0.  Box  1201 
San  Jose,  CA  95108 

Mr.  J.  Appelblatt  (01) 

Director  of  Engineering 
Cadillac  Gauge  Company 
P.0.  Box  1027 
Warren,  MI  48090 

Chrysler  Corporation  (02) 

Mobility  Research  Laboratory, 
Defense  Engineering 
Department  6100 
P.0.  Box  751  ' 

Detroit,  MI  48231 


CALSPAN  Corporation 
Box  235 
Library 

4455  Benesse  Street 
Buffalo,  NY  14221 

SEM,  (01) 

F  orsvaretsf  orskningsanstal t 
Avd  2 

Stockholm  80,  Sweden 

Mr.  Hedwig  (02) 

RU  III/6 

Ministry  of  Defense 
5300  Bonn,  Germany 

Foreign  Science  &  Tech  (01) 
Center 

220  7th  Street  North  East 
ATTN:  AMXST-GEI 

Mr.  Tim  Nix 

Charlottesville,  VA  22901 

General  Research  Corp  (01) 
7655  Old  Springhouse  Road 
Westgate  Research  Park 
ATTN:  Mr.  A.  Viilu 

McLean,  VA  22101 

Commander  ( 01 ) 

US  Army  Developmant  and 
Readiness  Command 
5001  Eisenhower  Avenue 
ATTN:  Dr.  R.S.  Wiseman 
Alexandria,  VA  22333 


President  (02) 

Army  Armor  and  Engineer  Board 
Fort  Knox,  KY  40121 

Commander  (01) 

US  Army  Arctic  Test  Center 
APO  409 

Seattle,  WA  98733 

Commander  (02) 

US  Army  Test  &  Evaluation 
Command 

ATTN:  AMSTE-BB  and  AMSTE-TA 
Aberdeen  Proving  Ground,  MD 
21005 

Commander  (01) 

US  Army  Armament  Research 
and  Development  Command 
ATTN :  Mr .  Rubin 

Dover,  NJ  07801 

Commander  (01) 

US  Army  Yuma  Proving  Ground 
ATTN :  STEYP-RPT 

Yuma,  AZ  85364 

Commander  (01) 

US  Army  Natic  Laboratories 
ATTN:  Technical  Library 

Natick,  MA  01760 

Director  (01) 

US  Army  Human  Engineering  Lab 
ATTN:  Mr.  Eckels 
Aberdeen  Proving  Ground,  MD 
21005 

Director  (02) 

US  Army  Ballistic  Research  Lab 
Aberdeen  Proving  Ground,  MD 
21005 

Director  (02) 

US  Army  Materiel  Systems 

Analysis  Agency 

ATTN:  AMXSY-CM 

Aberdeen  Proving  Ground,  MD 

21005 


Director  (02) 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  VA  22314 

US  Marine  Corps  (01) 

Mobility  &  Logistics  Division 
Development  and  Ed  Command 
ATTN:  Mr.  Hickson 
Quantico,  VA  22134 

Keweenaw  Field  Station  (01) 

Keweenaw  Research  Center 

Rural  Route  1 

P.0.  Box  94-D 

ATTN:  Dr.  Sung  M.  Lee 

Calumet,  MI  49913 

Naval  Ship  Research  &  (02) 

Dev  Center 

Aviation  &  Surface  Effects  Dept 
Code  161 

Washington,  D.C.  20034 

Director  (01) 

National  Tillage  Machinery  Lab 
Box  792 

Auburn,  AL  36830 

Director  (02) 

USDA  Forest  Service  Equipment 
Development  Center 
444  East  Bonita  Avenue 
San  Dimes,  CA  91773 

Engineering  Societies  (01) 

Library 

345  East  47th  Street 
New  York,  NY  10017 

Dr.  I.R.  Erlich  (01) 

Dean  for  Research 
Stevens  Institute  of  Technology 
Castle  Point  Station 
Hoboken,  NJ  07030 


